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ABSTRACT. The event of large losses plays an important role in credit risk. As these large losses are 
typically rare, and portfolios usually consist of a large number of positions, large deviation theory is 
the natural tool to analyze the tail asymptotics of the probabilities involved. We first derive a sample- 
path large deviation principle (LDP) for the portfolio's loss process, which enables the computation 
of the logarithmic decay rate of the probabilities of interest. In addition, we derive exact asymptotic 
results for a number of specific rare-event probabilities, such as the probability of the loss process 
exceeding some given function. 



1. Introduction 

For financial institutions, such as banks and insurance companies, it is of crucial importance to 
accurately assess the risk of their portfolios. These portfolios typically consist of a large number 
of obligors, such as mortgages, loans or insurance policies, and therefore it is computationally 
infeasible to treat each individual object in the portfolio separately. As a result, attention has 
shifted to measures that characterize the risk of the portfolio as a whole, see e.g. [DS03J for gen- 
eral principles concerning managing credit risk. The best known metric is the so-called value at 
risk, see |Jor06[, measuring the minimum amount of money that can be lost with a percent cer- 
tainty over some given period. Several other measures have been proposed, such as economic 
capital, the risk adjusted return on capital (RAROC), or expected shortfall, which is a coherent 
risk measure [ADEH99J. Each of these measure are applicable to market risk as well as credit 
risk. Measures such as loss given default (LCD) and exposure at default (EAD) are measures that 
purely apply to credit risk. These and other measures are discussed in detail in e.g. |MFE05|. 

The currently existing methods mainly focus on the distribution of the portfolio loss up to a given 
point in time (for instance one year into the future). It can be argued, however, that in many 
situations it makes more sense to use probabilities that involve the (cumulative) loss process, say 
{L(t) : t > 0}. Highly relevant, for instance, is the event that L(-) ever exceeds a given function 
£(■) (within a certain time window, for instance between now and one year ahead), i.e., an event 
of the type 

(1.1) {3t < T : L(t) > C(t)}. 

It is clear that measures of the latter type are intrinsically harder to analyze, as it does not suffice 
anymore to have knowledge of the marginal distribution of the loss process at a given point in 
time; for instance the event (Q) actually corresponds to the union of events {L(t) > £(t)}, for 
t < T, and its probability will depend on the law of L(-) as a process on [0, T]. 

In line with the remarks we made above, earlier papers on applications of large-deviation theory 
to credit risk, mainly address the (asymptotics of the) distribution of the loss process at a sin- 
gle point in time, see e.g. |DDD04| and [GKS07]. The former paper considers, in addition, also 
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the probability that the increments of the loss process exceed a certain level. Other approaches 
to quantifying the tail distribution of the losses have been taken by [LKSS01J, who use extreme- 
value theory (see [EKM97J for a background), |Gor02] and [MTB01 J, where the authors consider 
saddle point approximations to the tails of the loss distribution. Numerical and simulation tech- 
niques for credit risk can be found in e.g. [Gla04|. The first contribution of our work concerns 
a so-called sample-path large deviation principle (LDP) for the average cumulative losses for 
large portfolios. Loosely speaking, such an LDP means that, with L n (-) denoting the loss process 
when n obligors are involved, we can compute the logarithmic asymptotics (for n large) of the 
normalized loss process L n (-)/n being in a set of trajectories A: 



(1.2) lim -logP -£„(•) S A 

n^oo n \n 

we could for instance pick a set A that corresponds to the event \\.\) . Most of the sample-path 
LDPs that have been developed so far involve stochastic processes with independent or nearly- 
independent increments, see for instance the results by Mogul'skii for random walks |Mog76], de 
Acosta for Levy processes |Aco94], and Chang [Cha95| for weakly-correlated processes; results 
for processes with a stronger correlation structure are restricted to special classes of processes, 
such as Gaussian processes, see e.g. |Aze80|. It is observed that our loss process is not covered 
by these results, and therefore new theory had to be developed. The proof of our LDP relies 
on 'classical' large-deviation results (such as Cramer's theorem, Sanov's theorem, Mogul'skii's 
theorem), but in addition the concept of epi-convergence |Kal86J is relied upon. 

Our second main result focuses specifically on the event {1A) of ever (before some time horizon 
T) exceeding a given barrier function £(•). Whereas we so far considered, inherently impre- 
cise, logarithmic asymptotics of the type displayed in (|1.2b , we can now compute so-called exact 
asymptotics: we identify an explicit function f(n) such that f(n)/p n — > 1 as n — ► oo, where p n 
is the probability of our interest. As is known from the literature, it is in general substantially 
harder to find exact asymptotics than logarithmic asymptotics. The proof uses the fact that, after 
discretizing time, the contribution of just a single time epoch dominates, in the sense that there is 
a t* such that 

(1.3) P L n (t*) > C(n) / Pn 1, with Pn := P fat : ± L„(t) > ((t) 
This t* can be interpreted as the most likely epoch of exceeding £(•). 

Turning back to the setting of credit risk, both of the results we present are derived in a setup 
where all obligors in the portfolio are i.i.d., in the sense that they behave independently and 
stochastically identically. A third contribution of our work concerns a discussion on how to 
extend our results to cases where the obligors are dependent (meaning that they, in the terminology 
of BDDD04I , react to the same 'macro-environmental' variable, conditional upon which they are 
independent again). We also treat the case of obligor-heterogeneity: we show how to extend the 
results to the situation of multiple classes of obligors. 

The paper is structured as follows. In Section|2]we introduce the loss process and we describe the 
scaling under which we work. We also recapitulate a couple of relevant large-deviation results. 
Our first main result, the sample-path LDP for the cumulative loss process, is stated and proved 
in Section [3] Special attention is paid to, easily-checkable, sufficient conditions under which this 
result holds. As argued above, the LDP is a generally applicable result, as it yields an expression 
for the decay rate of any probability that depends on the entire sample path. Then, in Section|4] 
we derive the exact asymptotic behavior of the probability that, at some point in time, the loss 
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exceeds a certain threshold, i.e., the asymptotics of p n , as defined in 1 11. 3b . After this we derive a 
similar result for the increments of the loss process. Eventually, in Section^ we discuss a number 
of possible extensions to the results we have presented. Special attention is given to allowing 
dependence between obligors, and to different classes of obligors each having its own specific 
distributional properties. 

2. Notation and Definitions 

The portfolios of banks and insurance companies are typically very large: they may consist of 
several thousands of assets. It is therefore computationally impossible to estimate the risks for 
each element, or obligor, in a portfolio. This explains why one attempts to assess the aggregated 
losses resulting from defaults, e.g. bankruptcies, failure to repay loans or insurance claims, for 
the portfolio as a whole. The risk in the portfolio is then measured through this (aggregate) loss 
process. In the following sections we introduce the loss process and the portfolio constituents 
more formally. 

2.1. Loss Process. Let (fi, T , P) be the probability space on which all random variables below are 
defined. We assume that the portfolio consists of n obligors and we denote the default time of 
obligor i by t%. Further we write Ui for the loss incurred on a default of obligor i. We then define 
the cumulative loss process L n as 

n 

(2.1) L n (t) :=Y^UiZi{t), 

i=l 

where Zi(t) — 1 { T . < t j is the default indicator of obligor i. We assume that the loss amounts Ui > 
are i.i.d., and that the default times r, > are i.i.d. as well. In addition we assume that the loss 
amounts and the default times are mutually independent. In the remainder of this paper U and 
Z(t) denote generic random variables with the same distribution as the Ui and Zi(t), respectively. 

Throughout this paper we assume that the defaults only occur on the time grid N; in Section [5] 
we discuss how to deal with the default epochs taking continuous values. In some cases we 
explicitly consider a finite time grid, say {1,2,..., N}. The extension of the results we derive to 
a more general grid {0 < t± <ti < ... < ijv} is completely trivial. The distribution of the default 
times, for each j, is denoted by 

(2.2) Pr .= P(r = j), 

(2.3) F r .= F(T<j)=Y,Pi- 

i=l 

Given the distribution of the loss amounts Ui and the default times Tj, our goal is to investigate the 
loss process. Many of the techniques that have been developed so far, first fix a time T (typically 
one year), and then stochastic properties of the cumulative loss at time T, i.e., L n (T), are studied. 
Measures such as value at risk and economic capital are examples of these 'one-dimensional' 
characteristics. Many interesting measures, however, involve properties of the entire path of the 
loss process rather than those of just one time epoch, examples being the probability that that 
L„( ) exceeds some barrier function £(•) for some t smaller than the horizon T, or the probability 
that (during a certain period) the loss always stays above a certain level. The event corresponding 
to the former probability might require the bank to attract more capital, or worse, it might lead 
to the bankruptcy of this bank. The event corresponding to the latter event might also lead to the 
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bankruptcy of the bank, as a long period of stress may have substantial negative implications. 
We conclude that having a handle on these probabilities is therefore a useful instrument when 
assessing the risks involved in the bank's portfolios. 

As mentioned above, the number of obligors n in a portfolio is typically very large, thus pro- 
hibiting analyses based on the specific properties of the individual obligors. Instead, it is more 
natural to study the asymptotical behavior of the loss process as n — > oo. One could rely on a 
central-limit-theorem based approach, but in this paper we focus on rare events, by using the 
theory of large deviations. 

In the following subsection we provide some background of large-deviation theory, and we de- 
fine a number of quantities that are used in the remainder of this paper. 

2.2. Large Deviation Principle. In this section we give a short introduction to the theory of large 
deviations. Here, in an abstract setting, the limiting behavior of a family of probability measures 
{/i rl } on the Borel sets B of a metric space (X, d) is studied, asrn oo. This behavior is referred 
to as the Large Deviation Principle (LDP), and it is characterized in terms of a rate function. The 
LDP states lower and upper exponential bounds for the value that the measures p n assign to sets 
in a topological space X. Below we state the definition of the rate function that has been taken 
from ITDZ981 . 

Definition 2.1. A rate function is a lower semicontinuous mapping I : X — > [0, oo], for all a G [0, oo) 
the level set ^i(a) :— {x\ I(x) < a} is a closed subset of X. A good rate function is a rate function for 
which all the level sets are compact subsets of X. 

With the definition of the rate function in mind we state the large deviation principle for the 
sequence of measure {p n }- 

Definition 2.2. We say that {p n } satisfies the large deviation principle with a rate function !(■) if 

(i) (Upper bound) for any closed set F C X 

(2.4) limsup — \ogp n (F) < — inf I(x). 

n-too n x£F 

(ii) (Lower bound) for any open set G C X 

(2.5) lim inf - log p n (G) > - inf I(x) . 

n^oo n i£G 

We say that a family of random variables X — {X n }, with values in X, satisfies an LDP with rate function 
Ix{) iff the laws {p*} satisfy an LDP with rate function Ix, where p* is the law of X n . 

The so-called Fenchel-Legendre transform plays an important role in expressions for the rate 
function. Let for an arbitrary random variable X, the logarithmic moment generating function, 
sometimes referred to as cumulant generating function, be given by 

(2.6) A x (6) := logM x (9) = logE [e ex ] < oo, 

for 9 £ R. The Fenchel-Legendre transform of Ax is then defined by 

(2.7) A* x (x) :=sup(fe-A x (0)). 

e 

We sometimes say that A^ is the Fenchel-Legendre transform of X. 
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The LDP from Definition |Z2] provides upper and lower bounds for the log-asymptotic behavior 
of measures \i n . In case of the loss process l|2.1b , fixed at some time t, we can easily establish an 
LDP by an application of Cramer's theorem (Theorem IA.lt . This theorem yields that the rate 
function is given by A^ Z (t) (")/ where ^jjz(t) (■) is the Fenchel-Legendre transform of the random 
variable UZ(t). 

The results we present in this paper involve either A^( ) (Section [3), which corresponds to i.i.d. 
loss amounts Ui only or A£, z ^(-) (Section[4), which corresponds to those loss amounts up to 
time t. In the following section we derive an LDP for the whole path of the loss process, which 
can be considered as an extension of Cramer's theorem. 



3. A Sample-Path Large Deviation Result 

In the previous section we have introduced the large deviation principle. In this section we derive 
a sample-path LDP for the cumulative loss process (|2.1|> . We consider the exponential decay of 
the probability that the path of the loss process L n (-) is in some set A, as the size n of the portfolio 
tends to infinity. 



3.1. Assumptions. In order to state a sample-path LDP we need to define the topology that we 
work on. To this end we define the space of all nonnegative and nondecreasing functions on 
T N = {1,2,...,N}, 

y := {/ : T N -> R+| < f t < f i+1 for * < N}. 

This set is identified with the space M< := {x £ K. N \ < Xi < Xi+i for i < N}. The topology on 
this space is the one induced by the supremum norm 

I l/i loo = . max \fi\. 

t—l,...,N 

As we work on a finite-dimensional space, the choice of the norm is not important, as any other 
norm on would result in the same topology. We use the supremum norm as this is convenient 
in some of the proofs in this section. 

We identify the space of all probability measures on Tjv with the simplex $: 



(3.1) $ := U e 



y~] ^ = 1, tfi > for i < N 



For a given (p £ $ we denote the cumulative distribution function by tp, i.e., 

i 

(3.2) ipi = ^2cpj, fori < N; 

3=1 

note that tp G 5? and ipN = L 

Furthermore, we consider the loss amounts Ui as introduced in Section IZll a^e$ with cdf tp, 
and a sequence of ip n £ $, each with cdf ip n , such that ip n — > cp as n oo, meaning that ipf — > (pi 



6 



V. LEIJDEKKER, M. MANDJES, AND P. SPREIJ 



for all i < N. We define two families of measures (fi n ) and (v n ), 

( , , , N 

(3.3) ,i„(A):=P 



(3.4) v n {A) := 




S A 



e A 



i=l 



where A e B := B(S, N ) and [x] := supjfc e N\k < x}. Below we state an assumption under 
which the main result in this section holds. This assumption refers to the definition of exponential 
equivalence, which can be found in Definition E21 



Assumption 3.1. Assume that tp n — > ip and moreover that the measures \i n and v n as defined in (|3.3 
and (|3.4[) , respectively, are exponentially equivalent. 



^From Assumption 13.11 we learn that the differences between the two measures /i„ and v n go 
to zero at a 'superexponential' rate. In the next section, in Lemma [3.41 we provide a sufficient 
condition, that is easy to check, under which this assumption holds. 



3.2. Main Result. The assumptions and definitions in the previous sections allow us to state 
the main result of this section. We show that the average loss process satisfies a large deviation 
principle as in l|2.2|> . This principle allows us to approximate a large variety of probabilities related 
to the average loss process, such as the probability that the loss process stays above a certain time 
dependent level or the probability that the loss process exceeds a certain level before some given 
point in time. 

Theorem 3.2. With $ as in (|3.3[) and under Assumption ^. 1\ the average loss process, L n (-)/n satisfies 
an LDP with rate function Iu iP . Here, for x £ R< , Iu, P is given by 

(3.5) fe.W-«E*(^S) 
with Axi :— Xi — x^\ and xq := 0. 



Observing the rate function for this sample path LDP, we see that the effects of the default times 
Ti and the loss amounts Ui are nicely decomposed into the two terms in the rate function, one 
involving the distribution of the default epoch r (the 'Sanov term', cf. [DZ98 Thm. 6.2.10]), the 
other one involving the incurred loss size U (the 'Cramer term', cf. [DZ98 Thm. 2.2.3]. Observe 
that we recover Cramer's theorem by considering a time grid consisting of a single time point, which 
means that Theorem 13.21 extends Cramer's result. We also remark that, informally speaking, the 
optimizing (p 6 $ in d3.5b can be interpreted as the 'most likely' distribution of the loss epoch, 
given that the path of L n (-)/n is close to x. 

As a sanity check we calculate the value of the rate function Iu iP (x) for the 'average path' of 
L n { )/n, given by x* — E [U] Fj for j < N, where Fj is the cumulative distribution of the default 
times as given in 1)2. 3b ; this path should give a rate function equal to 0. To see this we first remark 
that clearly Ijj jP (x) > for all x, since both the Sanov term and the Cramer term are non-negative. 
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This yields the following chain of inequalities: 



N 



< I UiP (x*) = mf ^ Vi [ log ( g ) + A* 



< EMMS 



E|C/]Pj 



= 53piA&(E[Cn) = A&(E[l7]) = 0, 

i=l 

where we have used that for E [U] < oo, it always holds that A^(E [£/]) = cf. [DZ98, Lemma 
2.2.5]. The inequalities above thus show that, if the 'average path' x* lies in the set of interest, 
then the corresponding decay rate is 0, meaning that the probability of interest decays subexpo- 
nentially. 

In the proof of Theorem 13.21 we use the following lemma, which is related to the concept of epi- 
convergence, extensively discussed in |Kal86|. 



I, with D c K m compact. Assume that for all x £ D and for all x n 



Lemma 3.3. Let f n , / : D 

in D we have 

(3.6) lim sup /„(#„) < f(x). 

n — >oo 

Then we have 

limsup sup f n (x) < sup f(x). 

Proof. Let /* = sup xeD f n (x), f* = sup xeD f(x). Consider a subsequence f* k -> limsup,^^ /*. 
Let e > and choose x Uk such that /* < f nk {%n k ) + £ for all k. By the compactness of D, there 
exists a limit point x such that along a subsequence x nk . — > x. By the hypothesis d3.6|> we then 
have 

lim sup /* < lim f nk . (x nk . ) + e = f(x) + e < f* + e. 

n— >oo 3 

Let e { to obtain the result. □ 

Proof of Theorem l3.2l We start by establishing an identity from which we show both bounds. 
We need to calculate the probability 



L n {-) e A = 



1 



1 



-L n (l),...,-L n (N)\ e A , 
n n J J 

for certain A E B. For each point j on the time grid we record by the 'default counter' 
K n .j € {0, . . . , n} the number of defaults at time j: 

K n ,j '■ = #{« S {1, . . . ,n} | Tj = j}. 

These counters allow us to rewrite the probability to 

K n ,x+...+K n , K 

E ^ 



iL n (.) e A ) = e 

n 




1 

n 



(3.7) 



E w-' 



fci for i < N) x 
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where mj := Yl]=i kj an d the loss amounts Uj have been ordered, such that the first U(j) corre- 
spond to the losses at time 1, etc. 

Upper bound. Starting from Equality 03.7b , let us first establish the upper bound of the LDP. To 
this end, let F be a closed set and consider the decay rate 

(3.8) limsup - logP ( -in(-) e F 

An application of Lemma |A3] together with Il3.7b , implies that 0.8b equals 
limsup - logP ( -in(-) e ^ 



(3.9) = limsup max — 



io g p[ = i < .v ! i 



n n 



V 



/ v iV 



G F 



Next we replace the dependence on n in the maximization by maximizing over the set $ as in 
03 . lb - In addition, we replace the ki in d3.9b by 



(3.10) n>i 



[nipi] - [nipi 



where the ipi have been defined in (|3.2[) . As a result, 03. 8b reads 



(3.11) limsup sup — 



log I 



K„ 



<pn,u i<N)+\o, 




Note that 03.8b equals 03.11b , since for each n and vector (k±, . . . , k^) e N N , with $^ i=1 fci = n, 
there is a </? 6 $ with Lp i = fci/rt. On the other hand, we only cover outcomes of this form by 
rounding off the (fii . 

We can bound the first term in this expression from above using Lemma fA.51 which implies that 
the decay rate (13.8b is majorized by 



lim sup sup 

n — >oo 



JV 



^2 <Pn,i log 



<Pr. 




S F 



Now note that calculating the limsup in the previous expression is not straightforward due to the 
supremum over $. The idea is therefore to interchange the supremum and the limsup, by using 
Lemma [33] To apply this lemma we first introduce 



m ■= 



i=i v Fi / y\ j=i 




and note that $ is a compact subset of R". We have to show that for any sequence (p n — > i/3 
(3.12) limsup f n (<p n ) </(¥>), 



SAMPLE-PATH LARGE DEVIATIONS IN CREDIT RISK 



such that the conditions of Lemma f3.3l are satisfied. We observe, with ipf as in (|3.2|l and 0" as in 
(|3.10b with ip replaced by <p n , that 



limsup/„(^") < limsup - ^ tp^ log 

n — >oo n — >oo \ . ^ \ Pi 



lim sup — log 

n — >oc Tl 




Since (p n — > and since <p™ 4 differs at most by 1/n from <p™, it immediately follows that ip™ i — > 99 ^ 
For an arbitrary continuous function g we thus have g (<£>™ J — > g(tpi). This implies that 



Vn,i \ \ V~"> , / fi 



(3.13) limsup - V n ,j log ( I = - V log ( — 
Inequality (|3.12l l is established once we have shown that 



(3.14) limsup ilogpf^E 1 ^)) ^J^-i^E^(^) 



By Assumption 13.11 we can exploit the exponential equivalence together with Theorem IA.7I to 
see that (13.141 holds as soon as we have that 

/ , , , x A 



lim sup — log I 

n — »oc Tl 



\ 



;'|>0>) €F <-m|>a & (^) 



But this inequality is a direct consequence of Lemma IA.6I and we conclude that Q.14|l holds. 
Combining d3.13|) with (|3.14l l yields 



limsup /„(<?") < - V^log ( — ) - inf ( — ) 

n-»oo , =1 2:6 i=i \ / 



so that indeed the conditions of Lemma [373] are satisfied, and therefore 

limsup sup f n (<p) < sup f(<p) = sup ( - inf V tpi (log ( — ] + AJj ( ] ) 

n^oo v e* v g* y x ^ F ~[ \ \Pi J \ fi J J J 

&M5> flog fe) + (%fj) = - inf Wx) 



This establishes the upper bound of the LDP. 

Lower bound. To complete the proof, we need to establish the corresponding lower bound. Let G 
be an open set and consider 

(3.15) lim inf - logP \-L n {-) e G 



We apply Equality (|3.7l l to this liminf , with A replaced by G, and we observe that this sum is larger 
than the largest term in the sum, which shows that (where we directly switch to the enlarged 
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space $) the decay rate ([3.151 1 majorizes 



1 



liminf sup - logP -if n ,i = 0n,i, i < AT +lo. 



V 



)N 



e G 



Observe that for any sequence of functions h n (-) it holds that liminf„ sup,,, h n (x) > liminf„ h n (x) 
for all x, so that we obtain the evident inequality 

lim inf sup h n (x) > sup liminf h n (x). 

n — >oo x x n — >oo 

This observation yields that the decay rate of interest I l3.15b is not smaller than 



(3.16) sup 

V6* 



1 



V 



1 



liminf — logP —K n i = <p n ,i, i < N \ + liminf — log 



1 




where we have used that lim inf n (x n + y n ) > lim inf „ x„ + lim inf n y n . We apply Lemma [A.5I to 
the first liminf in ([3.161 1, leading to 

liminf - logP ( -K n i = <p ni , i < N ) > liminf [ — V* <p n i log f i £lhl \ _ — \og(n + 1) ] 

n— >oo U \n ) n— >oo \ f— \ Pi / Tl J 

(3.17) =-E^loggj, 

since log(n + l)/rt^0asrt^oo. The second liminf in d3.16|l can be bounded from below by an 
application of Lemma |Aj6] Since G is an open set, this lemma yields 

/ , , , x A 



lim inf — log I 

n— >oo ft 



\ 



' M>i] \ 



eG > " inf E^ At 



Upon combining these two results, we see that we have established the lower bound: 



liminf -logP -in(-) S G 

n — »oo 77, 



> - inf inf V <^ ( log 



inf Iu,p{x)- 



This completes the proof of the theorem. 



□ 



In order to apply Theorem 13.21 one needs to check that Assumption 13.11 holds. In general this 
could be a quite cumbersome exercise. In Lemma 13.41 below we provide a sufficient, easy-to- 
check condition under which this assumption holds. 



Lemma 3.4. Assume that for all 9 eM : Au(9) < oo. Then Assumption \3.1\ holds. 

Remark 3.5. The assumption we make in Lemma [3741 i.e., that the logarithmic moment generating 
function is finite everywhere, is a common assumption in large deviations theory. We remark that 
for instance Mogul'skii's theorem [DZ98. Thm. 5.1.2] also relies on this assumption; this theorem 
is a sample-path LDP for 

[nt] 



Y n [t) :=-Vl„ 

11 — ^ 
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on the interval [0, 1]. In Mogul'skii's result the Xi are assumed to be i.i.d; in our model we have 
that L n (t) = Yn=i UiZi(t)/n, so that our sample-path result clearly does not fit into the setup of 
Mogul'skii's theorem. Q 



Remark 3.6. In Lemma 
condition is 



(3.18) i im MM. = oo. 



it was assumed that Ajj(9) < oo, for all 9 £ R, but an equivalent 



In other words: this alternative condition can be used instead of the condition stated in Lemma 
13.41 To see that both requirements are equivalent, make the following observations. In Lemma 
IA.4I it is shown that l|3.18b is implied by the assumption in Lemma 13.41 In order to prove the 
converse, assume that l|3.18t holds, and that there is a < #o < oo for which Ajj(0) = oo. Without 
loss of generality we can assume that Ajj(0) is finite for 9 < do and infinite for 9 > 0$. For 
x > E [U], the Fenchel-Legendre transform is then given by 

A* u {x)= sup (9x-Au{9)). 

o<e<e 

Since U > and A^O) = 0, we know that A v (9) > for < 9 < 9 Q , and hence 



AM*) 



< 



which contradicts the assumption that this ratio tends to infinity as x 
the equivalence. 



oo, and thus establishing 







Proof of Lemma [3741 Let tp n — » p for some sequence of tp n G $ and tp e $. We introduce two 
families of random variables {Y n } and {Z n }, 

I i l !!*i J \ 

Y n 

which have laws /i„ and v n , respectively, as in 03.3b — J3.4|> . Since (p n — > ip we know that for any 
e > there exists an M £ such that for all n > M £ we have that max^ \<p? — ipi\ < e/N, and thus 
-tpi\<e. 

We have to show that for any 5 > 

"\Y n - Z n \\oo > S) = -oo. 
For i < N consider the absolute difference between Y n , and Z n j, i.e., 



Iim sup — log I 

n — >oc ^ 



(3.19) \Y n<i -Z n 



[mpi] 



iipi | < tie, which yields the upper bound, 



i=i j'=i 
Next we have that for any n > M e it holds that \mp\ 

\[n^]-[niPi]\ < [ne]+2, 

since the rounded numbers differ at most by 1 from their real counterparts. This means that the 
difference of the two sums in d3.19b can be bounded by at most [ne] + 2 elements of the Uj, which 
are for convenience denoted by U* . Recalling that the Uj are nonnegative, we obtain 

-, [™d i [ne]+2 

— n ^3 ' 



3=1 



3 = 1 
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Next we bound the probability that the difference exceeds 6, by using the above inequality: 

/, M+2 \ 

P (\\Y n - Z n \\ > 5) < p(- £ U*>5\ < E[cxp(9U 1 )} [n£]+2 e- n5e , 



where the last inequality follows from the Chernoff bound [DZ98 . Eqn. (2.2.12)] for 9 > 0. Taking 
the log of this probability, dividing by n, and taking the limsup on both sides results in 

limsup-logPdl^-^IU >8) <eA u (9)-59. 

n — >oo Tl 

By the assumption, Aj/(0) < oo for all 9. Thus, e — > yields 

limsup -logP(||y„ - > S) < -59. 

n — >oo 

As 9 was arbitrary, the exponential equivalence follows by letting 6 — ► oo. □ 
We conclude this section with some examples. 

Example 3.7. Assume that the loss amounts have finite support, say on the interval [0, u]. Then 
we clearly have 

A u (9) = logE [e eu ] <9u<oo. 

So for any distribution with finite support, the assumption for Lemma |3~4| is satisfied, and thus 
Theorem l3.2l holds. Here, the i.i.d. default times, Tj, can have an arbitrary discrete distribution on 
the time grid {1, . . . , N}. 

In practical applications one (always) chooses a distribution with finite support for the loss 
amounts, since the exposure to every obligor is finite. Theorem 13.21 thus clearly holds for any 
(realistic) model of the loss given default. 

An explicit expression for the rate function l|3.5t , or even the Fenchel-Legendre transform, is usu- 
ally not available. On the other hand one can use numerical optimization techniques to calculate 
these quantities. 

We next present an example to which Lemma l3.4l applies. 

Example 3.8. Assume that the loss amount U is measured in a certain unit, and takes on the values 
u, 2u, . . . for some u > 0. Assume that it has a distribution of Poisson type with parameter A > 0, 
in the sense that for i = 1,2,... 



—A 



A' 



P (U = (i + l)u) = e ,. 

i\ 

It is then easy to check that Ajj (9) = 9u + A (e 0u — l), being finite for all 9. Further calculations 
yield 



^w-e-'Mje-oH:- 1 )^ 



for all x > u, and oo otherwise. Dividing this expression by x and letting x — > oo, we observe 
that the resulting ratio goes to oo. As a consequence, Remark 13.61 now entails that Theorem 13.21 
applies. It can also be argued that for any distribution U with tail behavior comparable to that of 
a Poisson distribution, Theorem [372] applies as well. 
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4. Exact Asymptotic Results 

In the previous section we have established a sample-path large deviation principle on a finite 
time grid; this LDP provides us with logarithmic asymptotics of the probability that the sample 
path of L n (-)/n is contained in a given set, say A. The results presented in this section are dif- 
ferent in several ways. In the first place, we derive exact asymptotics (rather than logarithmic 
asymptotics). In the second place, our time domain is not assumed to be finite; instead we con- 
sider all integer numbers, N. The price to be paid is that we restrict ourselves to special sets A, 
viz. those corresponding to the loss process (or the increment of the loss process) exceeding a 
given function. We work under the setup that we introduced in Section l2Tl 

4.1. Crossing a Barrier. In this section we consider the asymptotic behavior of the probability 
that the loss process at some point in time is above a time-dependent level £. More precisely, we 
consider the set 

(4.1) A := {/ : T -> R+| 3t e T : f(t) > ((t)} , 
for some function ((t) satisfying 

(4.2) C(t) > E [UZ(t)\ = E [U] F t for all t e T, 

with F t as in (|2.3b . If we would consider a function £ that does not satisfy (|4.2) l, we are not in a 
large-deviations setting, in the sense that the probability of the event {L n {-)/n <E A} converges to 
1 by the law of large numbers. In order to obtain a more interesting result, we thus limit ourselves 
to levels that satisfy (|4.2b . For such levels we state the first main result of this section. 

Theorem 4.1. Assume that 

(4.3) there is a unique t* e T such that Iuz{t*) = niin/^z^), 
and that 

(4.4) liminf > 0, 

t^oo logt 

where I uz (t) = su Pe {6Q(i) - A UZ (t)(0)} = A* z(t) (C(t)) . Then 

for A as in (|4.ip and a* is such that A.' uz , tt Ja*) = ((t*). The constant C* follows from the Bahadur-Rao 
theorem (Theorem \A.8h with C* — Cuz(t*)x(t*)- 

Before proving our result, which will rely on arguments similar to those in [LM99J, we first dis- 
cuss the meaning and implications of Theorem 14.11 In addition we reflect on the role played by 
the assumptions. We do so by a sequence of remarks. 

Remark 4.2. Comparing Theorem 14.11 to the Bahadur-Rao theorem (Theorem IA.8K we observe 
that the probability of a sample mean exceeding a rare value has the same type of decay as 
the probability of our interest (i.e., the probability that the normalized loss process L n (-)/n ever 
exceeds some function £). This decay looks like C e^ nI / ^fri, for positive constants C and /. This 
similarity can be explained as follows. 
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First observe that the probability of our interest is actually the probability of a union events. 
Evidently this probability is larger than the probability of any of the events in this union, and 
hence also larger than the largest among these: 



(4.6) P (-L n (-) € A) > supP (-L n (t) > C(i) N 



Theorem 14. II indicates that the inequality in d4.6b is actually tight (under the conditions stated). 
Informally this means that the contribution of the maximizing t in the right-hand side of J4.6b , 
say t* , dominates the contributions of the other time epochs as n grows large. This essentially 
says that given that the rare event under consideration occurs, with overwhelming probability it 
happens at time t* . <J) 



As is clear from the statement of Theorem 14.11 two assumptions are needed to prove the claim; 
we now briefly comment on the role played by these. 

Remark 4.3. Assumption d4.3b is needed to make sure that there is not a time epoch t, different 
from t* t having a contribution of the same order as t*. It can be verified from our proof that if the 
uniqueness assumption is not met, the probability under consideration remains asymptotically 
proportional to e~ nI /y/n, but we lack a clean expression for the proportionality constant. 

Assumption (|4.4b has to be imposed to make sure that the contribution of the 'upper tail', that is, 
time epochs t G {t* + 1, t* + 2, . . .}, can be neglected; more formally, we should have 

3* E {t* + l,t* + 2,...} : -L n (t) > C(t)] =o(p (-L n {-) € A 

n J \ \n 

In order to achieve this, the probability that the normalized loss process exceeds C for large t 
should be sufficiently small. 



Remark 4.4. We now comment on what Assumption d4.4b means. Clearly, 

AtHj(t) (0) = logP (r < t) E [e eu ] + P (r > t) < logE [e eu ] , 

as t grows, 9 > 0; the limiting value is actually logE [e eu ] if t is non-defective. This entails that 

l uz (t) = A* z(tt) (C(t)) > AUC(t)) = sup (9at) - logE [e eu ]) . 

e 

We observe that Assumption l|4.4|l is fulfilled if liminft^oo A^j(((t))/ \ogt > 0, which turns out 
to be valid under extremely mild conditions. Indeed, relying on Lemma IA.41 we have that 
in great generality it holds that K\j{x)/x — > oo as x — > oo. Then clearly any ((t), for which 
lim inf t / log t > 0, satisfies Assumption @3j / since 

l imi nf^M =liminf AlipCW. 

t^oo lOgt t^co Q(t) lOgt 

Alternatively, if U is chosen distributed exponentially with mean A (which does not satisfy the 
conditions of Lemma |A. 4b , then A^(£) = At — 1 — log(At), such that we have that 

liminfM^ = A>0. 

t^oo lOgt 

Barrier functions C that grow at a rate slower than log t, such as log log t, are in this setting clearly 
not allowed. 
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Proof of Theorem l4.ll We start by rewriting the probability of interest as 



-L n (-)eA \ = ¥[3teT 



L n (t) 



>C(*) 



For an arbitrary point k in T we have 



(4.7) 



3t G T 



L n (t) 



> C(t) < P 3t < fc 



L n (t) 



>C(t) 



3t > k : > C(t) ) ■ 



We first focus on the second part in 04.7b . We can bound this by 



3t>k: ^>m)< E 



i=fc+l 



>C(» 



< V infE 



i=k+l 



exp \0Y,U S Z s (j) 
3=1 



where the second inequality is due to the Chernoff bound |DZ98, Eqn. (2.2.12)]. The indepen- 
dence between the Ui and Zi{t), together with the assumption that the Ui are i.i.d. and the Zi(t) 
are i.i.d., yields 



V infE 

^— ' 6>0 



i=k+l 



exp le^UjZjii) 



e -nme = g mf o nE[exp(^^(i))], 



<(i)8 



»=fc+l 3=1 



V exp [ -nsup(C(«)f- A [7Z(i) (6>) ] = V exp (-raI(/;j(C(0)) • 



i = k+l 



By (|4.4) l we have that 



t^oo iOgt 

for some (3 > (possibly oo). Hence there exists an m such that for alii > m 
(4.8) ii^(i)>alogi>Itrz(**)> 

where a = [3/2 (in case /3 = oo, any < a < oo suffices) and i* defined in 04 .3b . Choosing fc = m, 
we obtain by using the first inequality in 04.8b for n > 1/a 

x v i 

> exp(— nluz(((i))) < / exp (— nalogi) < exp ((— na + 1) logm) , 

* — ' z — ' na — 1 

i— m+1 z— m+1 

where the last inequality trivially follows by bounding the summation (from above) by an ap- 
propriate integral. Next we multiply and divide this by P (L n (t*)/n > C(O) anc ^ we a PPty the 
Bahadur-Rao theorem, which results in 

_J_ e (-tia+l)logm _ _J_ e (-na+l)logm P (L n (t*)/n > Cft*)) 

na -I nat-1 P (L n (t*)/n > £(**)) 



= p ( -L n (t*) > an 



na — 1 



1 - O [ - 

n 



,—n(a logm-Iu z (t*)) 



The second inequality in 04,8b yields a logm — Ijjz(t*) > 5, for some 5 > 0. Applying this 
inequality, we see that this bounds the second term in 04. 7b , in the sense that as n — > oo, 



3t> k 



^ > at)) I p ( -un > an 



lh 
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To complete the proof we need to bound the first term of fl4.7[l . where we use that k = m. For this 
we again use the Bahadur-Rao theorem. Next to this theorem we use the uniqueness of t* , which 
implies that for i < m and i ^t* there exists an e* > 0, such that 

Iuz(t*)+s* <I uz (i). 

This observation yields, with <r, such that A-' UZ U) ( ai ) = 

'L n (i) 



3t<m: 



< 



< 



L n (t*) > C(t*) J ll + O 
L n (i*) > C(**) ) fl + O 



L n (i*) > C(**) 1 + 



E 



— ^!7Z(i),C(i) 



1 + m x max 

z— l,...,m 



-nluz(ti) 



c* 



c, 



uz(i),c(i) 



(--)) 



Combining the above findings, we observe 

'£«(**) 



3t g T : > C(t) ) < 



Together with the trivial bound 



3t e T : > C(t) ] : 



1 + 



> C(**) 1 + 



this yields 



3teT: ^±>((t) 



> an 



>C(0 i + o - 



Applying the Bahadur-Rao theorem to the right-hand side of the previous display yields the 
desired result. □ 



4.2. Large Increments of the Loss Process. In the previous section we identified the asymptotic 
behavior of the probability that at some point in time the normalized loss process L„(-)/ n exceeds 
a certain level. We can carry out a similar procedure to obtain insight in the large deviations of the 
increments of the loss process. Here we consider times where the increment of the loss between 
time s and t exceeds a threshold £(s, t). More precisely, we consider the event 

(4.9) A:={f:TxT^R+\3s<teT: f(s,t)>£(s,t)}. 

Being able to deal with events of this type, we can for instance analyze the likelihood of the 
occurrence of a large loss during a short period; we remark that with the event 04. 1 b from the 
previous subsection, one cannot distinguish the cases where the loss is zero for all times before 
t and x > C(t) at time t, and the case where the loss is just below the level ( for all times before 
time t and then ends up at x at time t. Clearly, events of the d4.9b make it possible to distinguish 
between such paths. 

In order to avoid trivial results, we impose a condition similar to l|4.2|l . namely 

(4.10) £(8 ) t)>E[U] (F t — F s ) , 

for all s < t. The law of large numbers entails that for functions £ that do not satisfy this condition, 
the probability under consideration does not correspond to a rare event. 
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A similar probability has been considered in [DDD04J, where the authors derive the logarithmic 
asymptotic behavior of the probability that the increment of the loss, for some s < t, in a bounded 
interval exceeds a thresholds that depends only on t — s. In contrast, our approach uses a more 
flexible threshold, which depends on both times s and t, and in addition we derive the exact 
asymptotic behavior of this probability. 

Theorem 4.5. Assume that 

(4.11) there is a unique s* <i*eT such that Iuz{s* , £*) = mmluz(s,t), 

s<t 

and that 

(4.12) inf lim inf kl^A > 
s£T t^oo logt 

where I uz (s,t) = sup e (0£(s,f) - ^u{z{t)-z{ s )){0)) = A y(z (f )-z( s ))(£(M))- Then 

( 1 \ e -nIuz{s*,t*) C * / /j 

(4.13) P -L„(.) eA = = 1 + O 



for A as in ( 14.1 OP and a* is such that ^■' u ^ Z (t*)-z(s*))( a *) = £( s *> ^ e constant C* follows from the 
Bahadur-Rao theorem (Theorem \A.8h with C* = Cijrz(t*)-z(s*)), £(»*,{*)• 

Remark 4.6. A first glance at Theorem l4.5l tells us the obtained result is very similar to the result of 
Theorem 14. II The second condition, i.e., Inequality ( |4.12> , however, seems to be more restrictive 
than the corresponding condition, i.e., Inequality l|4.4) l, due to the infimum over s. This assump- 
tion has to make sure that the 'upper tail' is negligible for any s. In the previous subsection we 
have seen that, under mild restrictions, the upper tail can be safely ignored when the barrier 
function grows at a rate of at least logt. We can extend this claim to our new setting of large 
increments, as follows. 

First note that 

inf lim inf > inf lim inf A &p''». 

s£T t— >oo logi sGT i— >oo log t 

Then consider thresholds that, next to condition l|4.10t , satisfy that for all s 

(4.14) liminf^^>0. 

t^oo logt 

Then, under the conditions of Lemma fA.41 we have that 

(4.15) liminf ^t)) = A* p t)) tOrf = 

since the second term remains positive by l|4.14|l and the first term tends to infinity by Lemma 
IA.4I Having established i4.15|l for all s, it is clear that (|4.12|l is satisfied. 

The sufficient condition l|4.14b shows that the range of admissible barrier functions is quite sub- 
stantial, and, importantly, imposing l|4.12b is not as restrictive as it seems at first glance. 

Proof of Theorem l4.5l The proof of this theorem is very similar to that of Theorem l4.ll Therefore 
we only sketch the proof here. 

As before, the probability of interest is split up into a 'front part' and 'tail part'. The tail part 
can be bounded using Assumption < |4.12> ; this is done analogously to the way Assumption l|4.4l l 
was used in the proof of Theorem 14.11 The uniqueness assumption l|4.11|l then shows that the 
probability of interest is asymptotically equal to the probability that the increment between time 
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s* and t* exceeds £(s* , t*); this is an application of the Bahadur-Rao theorem. Another application 
of the Bahadur-Rao theorem to the probability that the increment between time s* and t* exceeds 
£(s* , t*) yields the result. □ 

5. Discussion and Concluding Remarks 

In this paper we have established a number of results with respect to the asymptotic behavior of 
the distribution of the loss process. In this section we discuss some of the assumptions in more 
detail and we consider extensions of the results that we have derived. 

5.1. Extensions of the Sample-Path LDP. The first part of our work, Section|3j was devoted to 
establishing a sample-path large deviation principle on a finite time grid. Here we modelled 
the loss process as the sum of i.i.d. loss amounts multiplied by i.i.d. default indicators. ^From a 
practical point of view one can argue that the assumptions underlying our model are not always 
realistic. In particular, the random properties of the obligors cannot always be assumed indepen- 
dent. In addition, the assumption that all obligors behave in an i.i.d. fashion will not necessarily 
hold in practice. Both shortcomings can be dealt with, however, by adapting the model slightly. 

A common way to introduce dependence, taken from [DDD04J, is by supposing that there is a 
'macro-environmental' variable Y to which all obligors react, but conditional on which the loss 
epochs and loss amounts are independent. First observe that our results are then valid for any 
specific realization y of Y. Denoting the exponential decay rate by r y , i.e., 



lim -logP ( -£„(•) E A 



+oo n \ n 



Y = y 



the unconditional decay rate is just the maximum over the r y ; this is trivial to prove if Y can attain 
values in a finite set only. A detailed treatment of this is beyond the scope of this paper. 

The assumption that all obligors have the same distribution can be relaxed to the case where we 
assume that there are m different classes of obligors (for instance determined by their defaultrat- 
ings). We further assume that each class i makes up a fraction <n of the entire portfolio. Then we 
can extend the LDP of Theorem l3.2l to a more general one, by splitting up the loss process into to 
loss processes, each corresponding to a class. Conditioning on the realizations of these processes, 
we can derive the following rate function: 

m N 

(5.1) Iu, P , m (x) := inf 'mf} } ai(f^ 

3=1 i=l 




where V x = { v e K+ x ^ | Yl?=i v l = Ax i for a11 i < N }> and $ ™ is the Cartesian product $ 



x . 



$ (to times), with $ as in l|3.1b . The optimization over the set V x follows directly from conditioning 
on the realizations of the per-class loss processes. We leave out the formal derivation of this 
result; this multi-class case is notationally considerably more involved than the single-class case, 
but essentially all steps carry over. 

In our sample-path LDP we assumed that defaults can only occur on a finite grid. While this 
assumption is justifiable from a practical point of view, an interesting mathematical question is 
whether it can be relaxed. In self-evident notation, one would expect that the rate function 



Iu,p,oo(x) := inf ifj ( log ( — ) + ^ 
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It can be checked, however, that the argumentation used in the proof of Theorem 13.21 does not 
work; in particular, the choice of a suitable topology plays an important role. 

If losses can occur on a continuous entire interval, i.e., [0, N], we expect, for a nondecreasing and 
differentiable path x, the rate function 

(5.2) I UM ( X ) := V (t) (log + A h («)) d*. 

where M is the space of all densities on [0, N] and p the density of the default time r. One can 
easily guess the validity of (|5.2b from d3.5b by using Riemann sums to approximate the integral. 
A formal proof, however, requires techniques that are essentially different from the ones used to 
establish Theorem l3.21 and therefore we leave this for future research. 

5.2. Extensions of the Exact Asymptotics. In the second part of the paper, i.e., Section |H we 
have derived the exact asymptotic behavior for two special events. First we showed that, under 
certain conditions, the probability that the loss process exceeds a certain time-dependent level, 
is asymptotically equal to the probability that the process exceeds this level at the 'most likely' 
time t*. The exact asymptotics of this probability are obtained by applying the Bahadur-Rao 
theorem. A similar result has been obtained for an event related to the increment of the loss 
process. One could think of refining the logarithmic asymptotics, as developed in Section 
to exact asymptotics. Note, however, that this is far from straightforward, as for general sets 
these asymptotics do not necessarily coincide with those of a univariate random variable, cf. 
IIMMNvU06l . 
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Appendix A. Background Results 

In this section we state a number of definitions and results, taken from [DZ98J, which are used in 
the proofs in this paper. 

Theorem A.l (Cramer). Let X; be i.i.d. real valued random variables with all exponential moments finite 
and let \x n be the law of the average S n — 2~2i=i Xi/n. Then the sequence {/in} satisfies an LDP with rate 
function A*(-), where A* is the Fenchel-Legendre transform of the X t . 



Proof. See for example 1DZ98J Thm. 2.2.3]. □ 

Definition A.2. We say that two families of measures {p n } and {u n } on a metric space (X, d) are expo- 
nentially equivalent if there exist two families of X -valued random variables {Y n } and {Z n } with marginal 
distributions {p. n } and {v n }, respectively, such that for all S > 

limsup-logP(d(A„,r„) > 8) = -oo. 

n — >oc n 

Lemma A.3. For every triangular array a l n > 0, n > 1, 1 < i < n, 

1 " 1 
lim sup — log a l n = lim sup max — log a l n . 

Proof. Elementary, but also a direct consequence of [DZ98 Lemma 1.2.15]. □ 
Lemma A.4. Let A(0) < oo for all 9 G K, then 

lim — j — j — = oo. 

\x\^oo \x\ 

Proof. This result is a part of 1PZ98J Lemma 2.2.20]. □ 

Lemma A.5. Let K n ,i be defined as K n j := ff{i e {1, . . . , n} \ n — j}. Then for any vector k € N N , 
such that J2iLi hi = n > we nave that 

[n + 1)- N exp (-nH(k \ p)) < P (K n = k) < exp (-nH(k \ p)) , 

where 

N k f k 

H{k\ P ) = Y,-M— 

\ n Pi 



and pi as defined in $2.2 



Proof. See [DZ98, Lemma 2.1.9]. 



□ 
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Lemma A.6. Define 

[nt] 

Z n (t) :=-Vjfi, 0<t<l, 
n 

i—l 

for an i.i.d. sequence ofW 1 valued random variables X t . Let //„ denote the law of Z n {-) in £00 ([0, 1])- For 
any discretization J = {0 < t\ < . . . < t\j\ < 1} and any f : [0, 1] — > R d , let pj(f) denote the vector 
(/(it))|=i € (K rf )l J L T/zen sequence of laws {[i n o pj 1 } satisfies the LDP in (M d )l J l roff/i the good 
rate function 

i J (z)=j2( t i- t i-i) A ' 

i=l 

where A* fs ffte Fenchel-Legendre transform of X\. 




Proof. See [DZ98, Lemma 5.1.8]. This lemma is one of the key steps in proving Mogul'skii's 
theorem, which provides a sample-path LDP for Z n {-) on a bounded interval. □ 

Theorem A.7. If an LDP with a good rate function /(•) holds for the probability measures {^„}, which 
are exponentially equivalent to {v n }, then the same LDP holds for {v n }. 



Proof. See JDZ9S Thm. 4.2.13]. □ 
Theorem A.8 (Bahadur-Rao). Let Xi be an i.i.d. real-valued sequence random variables. Then we have 



n 1 — ' / \/n \ \n 

v 2=1 / 



The constant Cx, q depends on the type of distribution of X\, as specified by the following two cases. 

(i) The law of ' X\ is lattice, i.e. for some xq, d, the random variable (X\ — xo)/d is (a.s.) an integer 
number, and d is the largest number with this property. Under the additional condition < 
P (Xi = q) < 1, the constant Cx, q is given by 

r d 



(1-e— V 2 ^ A xM 

where a satisfies K' x (a) = q. 
(ii) If the law of X\ is non-lattice, the constant Cx, q is given by 

Cx, q 



<V2tA£(<t)' 
with a as in case (i). 



Proof. We refer to HBR60II or IDZ98I Thm. 3.7.4] for the proof of this result. □ 
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